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ABSTRACT 

Recent X-ray observations have shown broad, skewed iron hne emission from Seyfert 1 
galaxies which is explained by the emission being fluorescence on a disk close to a black 
hole. During one interval, the line in MCG-6-30-15 was so broad and redshifted that 
a Kerr black hole is implied. We are therefore studying the effects of the Kerr metric 
on the line profile, and extending the work by Laor and Kojima which dealt only with 
extreme values of the spin parameter. Here we report that the spin parameter of the 
black hole in MCG-6-30-15 is high [a/M > 0.94), and invert the line profile to obtain 
the disk emissivity profile, which approximates a power-law. Continuum radiation 
returning to the disk because of the Kerr metric does not enhance the equivalent 
width of the line seen above 3 keV by more than about 20 per cent if the continuum 
source corotates with the disk. 
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1 INTRODUCTION 



X-ray observations of Seyfert 1 galaxies made with ASCA 
have shown that the iron emission lines discovered with 
GINGA (Pounds et al. 1990; Matsuoka et al. 1990) are broad 
(Fabian et al. 1994; Mushotzky et al. 1995; Reynolds 1997; 
Nandra et al. 1997). In particular, the line in MCG-6-30- 
15 was found to be both broad and skew from a 4.5 day 
long observation in 1994 (Tanaka et al. 1995). The line ex- 
tends about 2 keV below the rest frame energy of the line 
(6.4 keV) and only 0.3 keV above. The line profile is well 
fitted by that expected from fluorescence of matter on the 
surface of an accretion disk less than 20 gravitational radii 
(i.e. 20rg = 20G'M/c^) from a black hole of mass M (Tanaka 
et al. 1995; Fabian et al. 1995). Much of the skewness of the 
line is explained by gravitational redshift, and is the first 
clear evidence for the effects of strong gravity. 

Recent work examining variability in the line profile 
during the 4.5 day observation has revealed that it broad- 
ened still further during a deep minimum in the light curve 
(Iwasawa et al. 1996). The red (i.e. low energy) wing of the 
line extended further to the red and the blue (i.e. high en- 
ergy) wing disappeared. This behaviour can be explained if 
much of the emission originates from within 6rg (Iwasawa et 
al. 1996) . Since a disk around a non-spinning, Schwarzschild 



black hole does not extend within this region, it was con- 
cluded that the black hole must be spinning rapidly. In this 
case the Kerr metric applies and frame dragging causes the 
disk to extend inward of 6rg. 

Light bending effects, in particular the return of some 
disk radiation to the disk itself, can then become signifi- 
cant (see e.g. Cunningham 1976). In Iwasawa et al. (1996) 
it was proposed (without detailed calculation) that these ef- 
fects could enhance the strength (equivalent width) of the 
fluorescent line relative to to the continuum, emitted above 
the disk, as required by the observations. 

In this Letter we examine the line emission from a disk 
around a Kerr black hole, from the point of view of pro- 
file and of equivalent width, for the parameters relevant to 
MCG-6-30-15. Previously only the line profile of either a 
Schwarzschild black hole (Fabian et al. 1989) or a maxi- 
mal Kerr black hole (Laor 1991; Kojima 1991) have been 
available for fitting. Here we explicitly fit for the angular 
momentum parameter of the black hole, as well as the incli- 
nation angle of the normal to the disk to the line of sight. 
In addition, results of a free-form fit to the radial emissivity 
profile of the disk in the fiuorescent line are given, again for 
the first time. 
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2 DISKLINE MODEL USING THE KERR 
METRIC 

The propagation of radiation from a disk around a Kerr 
black hole has already been studied (e.g. Cunningham 1975, 
1976) and applied to model the observed profile of a line 
emission from the disk (Asaoka 1989; Kojima 1991; Laor 
1991; Bromley, Chen & Miller 1997). The main parameters 
to which the spectra are sensitive are the inclination angle 
of the observer from the axis of symmetry i, the radial emis- 
sivity profile of the disk e{r) and the angular momentum 
parameter of the hole a, = a/M, where a is the angular 
momentum and M is the black hole mass. (We note the 
model cannot discriminate a and M separately). We shall 
here investigate the effect of changes in i and a, , and show a 
derived emissivity profile across the disk, at fixed i and o*. 
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Figure 1. Predicted spectral line shapes for a Schwarzschild black 
hole (a, = 0). A is for i = 85°, B for i = 30°, and C for i = 0°. 



2.1 Assumptions and method of calculation 

• The black hole belongs to the Kerr family of solutions. 
The metric used is in the Boyer-Lindquist form. The grav- 
ity of the disk itself (and of the clouds responsible for the 
fluorescence) is negligible. 

• The accretion disk is assumed to be geometrically thin 
(i.e. thickness <^ smallest radius) and axially symmetric 
around the hole axis of rotation. 

• The matter responsible for the fluorescent line is sup- 
posed to lie on the surface of the disk and rotates on equato- 
rial circular geodesies. Emission starts at radius ri„, which 
is fixed to the radius of marginal stability Vms (Novikov & 
Thorne 1973), up to an arbitrary outer radius Vout- No signif- 
icant emission is expected far away from the centre since the 
line emissivity, described below, decreases rapidly (0.4% and 
0.002% of its maximum for radii of 15rg and lOOrg respec- 
tively if a, = 0.998). Previous studies (Tanaka et al. 1995) 
suggest an outer radius of ~ 10rg-16rg, so that we decided 
to fix this parameter to rout ~ 15rg. This constraint might 
have to be relaxed in future studies. 

• Except for the case whe re we explicitly fit for the emis- 
sivity profile (see Section 3.1 below), then in the rest frame of 
the emitting material (which we usually call the disk frame), 
the emissivity of the fluorescent component e(r) [defined as 
energy emitted per unit proper time per unit proper area] 
is assumed to follow the law described in equation (lib) 
of Page & Thorne (1974). This law is derived for the con- 
tinuum emission of the disk, but we suppose that the line 
flux (integrated across frequency in the emission rest frame) 
also follows this law. Since e(r) so defined vanishes at both 
r^ns and r = oo, we expect this emissivity profile to be more 
realistic than a simple power law model (Fabian et al. 1989). 

• We assume that the emitted fluorescent iron Ka line 
can be approximated by a 5-function in frequency, so that 
each emitted photon has an energy huo of 6.40 keV in the 
disk frame (Iwasawa et al. 1996). The emission is assumed 
to be locally isotropic. Thus, in the disk frame, the emitted 
intensity I^{r,v^) and flux Fe{r,i/e) are described by: 



Fe{r, Ve) = 7r/e(r, v^) = e{r)5{ve - ^o)- 



(1) 



• The fluorescent photons are free to reach the observer 
unless they cross the hole event horizon (lost radiation) or 
they return to the disk itself. (See Section ^ and Cunning- 
ham 1976). The observer is assumed to be located in the 




Figure 2. Same as Fig. ^ but for a» = 0.5. 



asymptotically flat region of space-time. In practice we use 
a distance of lOOOrg. 

The model then carries out a immerical integration of 
photon null geodesies backwards in time, from the observer 
to the disk, taking into account all the relativistic effects 
on both trajectories and redshift. The photons are emitted 
from the observer through a grid of equal solid angle pixels, 
which could be those of a CCD detector. At the end of the 
integration process, each image pixel is then associated with 
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Figure 3. Same as for Fig. but for a* = 0.998, corresponding 
to an extreme Kerr black hole (Thorne, 1974). 
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Figure 4. Predicted spectral line shapes at a fixed inclination 
i = 30°. D is for an extreme Kerr black hole (a, = 0.998), E for 
a* = 0.5, and F for a Schwarzschild black hole (a* = 0). 
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Figure 5. Contributions to the total lineshape from the receding 
(dashed) and approaching (solid) sections of the disk for i = 30° , 
a* = 0.998, Tin = lOrg and rout = 15rg. 



an emission radius re on the disk and a redshift z. We may 
compute the flux at the observer by using the phase space 
occupation number I^jv^ ^ which is invariant along the entire 
trajectory of the photon. Considering an interval v to v + 
Au in observed frequency, then the observed flux in a pixel 
corresponding to solid angle AQ at the observer, will be 



F{re 



AQ £(re) 
TvAu (1 + 2)''' 



(2) 



if the interval {u, u + Ai/) contains ^0/(1 + 2), and zero other- 
wise. Since redshift and flux are known for each pixel, images 
of those quantities can be plotted, such as those in Luminet 
(1979) or Bromley et al. (1997), but covering a wider range 
of parameters. Binning all f (re, z) into the redshift axis, en- 
ables one to compute the observed line proflle as seen by a 
distant observer. 
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Figure 6. Contours of probability versus i and a» are shown, 
calculated assuming a uniform prior on each parameter. 



2.2 The line profile as a function of spin 
parameter and inclination 

Figs. |l[ ^ and ^ show predicted spectral lines over a range 
of a, and i parameters and suggest that, for flxed emissivity 
profile, the overall line shape, in particular its width, is most 
sensitive to the inclination angle. For sufficiently high i (i.e. 
i > 25°) the spectrum is double-peaked. The low redshift 
component is mainly due to the receding part of the disk 
from the observer and the high redshift photons essentially 
come from the approaching half-part, as shown below. In the 
case of more face-on inclinations (i.e. i < 25°), the radial 
Doppler effect is not efficient enough and the observed line 
shape is a one-peak line with a weak blue component and a 
rather spread red contribution. Over the whole range of i, 
the red wing of the spectrum is being stretched by transverse 
Doppler and gravitational redshifts. This effect is boosted as 
the angular momentum a increases since the inner part of 
the disk then lies closer to the hole. Thus for a given i, 
the angular momentum controls the red wing extension (see 
Fig. For a given a*, i mainly constrains the upper limit 
of the blue part of the line (see Figs. |l[ |^ and 

The radial Doppler shifts are illustrated in Fig. |5|, where 
receding and approaching halves of the disk contributions 
are separated. An extreme black hole is used in Fig. |^ but 
with an emission area located between ri„ — lOrg and rout = 



15rg, in order to show a case where the two components are 
clearly separated in frequency. Such an appearance may be 
detectable if the X-ray emission is from a few localized flares 
which move around on the disk (see e.g. Iwasawa et al. 1996). 
Note the asymmetry of the two components, with the blue 
peak being much more intense than the red one. This blue 
boosting effect arises from the rather large redshift difference 
between the two halves. 



3 FITS TO THE LINE PROFILES OF 
MCG-6-30-15 

The evidence of strong gravity around a Kerr hole in MCG- 
6-30-15 has already been discussed in Iwasawa et al. (1996). 
This discussion included a comparison between the observed 
fluorescent Kc line at the period of lowest luminosity and 
a predicted Kerr lineshape corresponding to i = 30° and 
a, = 0.998 (Laor 1991). We report here an actual fit using 
our diskline model in order to obtain joint constraints on 
both a, and i. The data consist of 18 points lying between 
~ 3 keV and ~ 8 keV and are described in detail in Iwasawa 
et al. (1996). We have carried out a St assuming the 
emissivity profile of Page & Thorne as discussed above. The 
results are presented in Figs. ^ and ^. 

The contours give relatively strong evidence for an 
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Actual versus predicted line shape 



Figure 7. The line sliape corresponding to tlie best-fit values 
i = 29°^3'2i = 1—0 01 compared to the data points from 
Iwasawa et al. (1996). The error bounds on the best fit parameters 
are 68% confidence intervals obtained by marginalization with a 
uniform prior over the other parameter. 



inclination angle of ~ 25°-30° and strongly favour an ex- 
treme Kerr hole (a« > 0.94) rather than a Schwarzschild hole 
(a* = 0). The inclination i ~ 30° arises mainly because of 
the observed cut-off for energies higher than 6.5 keV, while 
the observed broad red wing necessitates the high value of 
a* . We note that a simple calculation shows that significant 
line emission below 3.8 keV requires that the disk extend 
within 1.9rg which in turn requires a, > 0.95. 



3.1 Fitting the emissivity profile 

The need for an a» value near 1 is so crucial in obtaining 
the highly redshifted emission seen in the data, that the 
observed line can strongly constrain the radial emissivity 
profile, e{r). This is because while a* near 1 allows a much 
smaller inner radius of the disk, it is only possible to take 
advantage of this if we weight the inner regions highly in 
forming the line profile. This weighting function is of course 
basically the emissivity profile. We have assumed above the 
Page & Thorne form, which (for the parameters of interest 
here) behaves approximately like r~^'* at large r and cuts 
off to zero at Vms- It is of interest to try to establish quan- 
titatively the e(r) profile favoured by the data itself, to see 
if this can reproduce something like a power law form, or 
even the cut off at rms, favoured theoretically. As a first 
test case, one might attempt to do this with the other im- 
portant parameters, the inclination angle and the black hole 
angular momentum parameter, held fixed. At a later stage 
one could then try to include fitting of all three of i, a* and 
a free-format e(r) simultaneously. We have written inver- 
sion software, incorporating maximum entropy regulariza- 
tion, which can achieve this, and this will be described fully, 
together with simulations, in a later paper. To show results 
of immediate interest in relation to MCG-6-30-15, we show 
in Figs. H and |^ the results obtained by this method with- 
out any regularization, using just positivity of the emission 
as the only constraint, and with i and a, fixed at values 
plausible on other grounds of 30° and 0.998 respectively. 
Such an approach is similar to the Lucy algorithm of Man- 
nuci, Salvati & Stanga (1992) in application to inverting the 
profiles of optical Balmer lines, except in the present case 
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Figure 8. Comparison of actual data points (1 sigma error bars 
shown) with the results in each bin predicted using the inferred 
emissivity profile (open circles). 



Inferred Disk Emissivity Profile 
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Figure 9. Log-log plot of inferred emissivity profile. 



the integral kernel is more complicated and has to be eval- 
uated numerically. Fig ^ shows the inferred emissivity pro- 
file as found in 20 radial bins, stretching from rms to 15rg 
in r. Some evidence for a power law form is evident, with 
a slope of perhaps ~ —3.5, although it obviously becomes 
noisy at the smaller flux levels. Fig. ^ compares the points 
predicted using this radial profile (and assuming i — 30° 
and a« = 0.998) with the observed data. It is evident that 
a degree of overfitting may be taking place, but that the 
results look sensible. 



4 THE ENHANCEMENT OF EQUIVALENT 
WIDTH BY RETURNING RADIATION 

Light bending effects become very powerful for highly rotat- 
ing holes and should contribute to enhance the line equiv- 
alent width. This bending effect is illustrated in Fig. |l^ by 
a polar diagram showing that for high inclination angles 
the emission strongly diverges from an isotropic case. The 
photons released from the inner part of the disk are prefer- 
entially bent toward the plane of rotation. 

Moreover, a substantial number of the geodesies re- 
turn to the disk itself — this is the returning radiation first 
treated by Cunningham (1976). To study the returning radi- 
ation, we used the same method of calculation as described 



in Section 2.1, but placing the observer at a given radius 
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Figure 10. Received power as a function of emitted radius (from 
left to right: r = 2rg, r = 8rg). The circular curve is for isotropic 
emission. 
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Figure 11. Returning radiation ratio : Fret(?'ret)/e(fi-et) (solid 
line). The dashed line shows the ratio of returning photons over 
the number of emitted photons; it is proportional to the enhance- 
ment in equivalent width (for a continuum power-law of energy 
index 1). 

Tret on the rotating disk. The result is shown in Fig. |ll] as 
the ratio of returning flux Fretij-ret) (coming from any part 
of the disk) over emitted flux e{rret) at the same point (as- 
suming again the Page & Thorne law). One can notice that 
for inner radii (i.e. Tret <~ l-4fg) this ratio is larger than 
unity. Two effects from the returning radiation should be 
involved in enlarging the equivalent- width of the fluorescent 
line relative to the primary power-law continuum: 

(i) More than half the continuum flux may hit the disk 
and produce fluorescent photons. 
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(ii) Some of the returning photons are strongly 
blueshifted (Fig. |l|). This effect is particularly enhanced 
on the inner parts of the disk. Then, since the power law 
slope of the continuum is negative (~ —2), at any given fre- 
quency, the continuum flux incident on the disk is higher 
than the continuum flux which reaches the observer. 

Since most of the line emission within r ~ 1.4rg emerges 
below 2 keV, the line enhancement due to returning radia- 
tion has little effect on the observed profile. An enhance- 
ment is instead needed at ~ 4 keV, or from r ~ 2rg, where 
returning radiation has little effect. Either the geometry of 
Martocchia & Matt (1996), in which the continuum source 
does not corotate with the disk, or some other effect (e.g. 
high ionization of iron in the disk) is needed for the observed 
increase in line strength (Iwasawa et al. 1996). 

More detailed calculations of line strength in both the 
geometry considered in this paper and that of Martocchia 
& Matt will be presented in a later paper (Lasenby et al. in 
preparation). This will cover a wider range of applications, 
including predicted disk images and some differences from 
previous results for line profiles from accretion disk will be 
highlighted. 
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Figure 12. The line profile as seen by an observer rotating with 
the disk at a radius of 1.5rg (solid line), or 2rg (dashed line). 
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